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ABSTRACT 

Two-dimensional  turbulent  v;akes  in  constant  adverse  and 
favorable  pressure  gradients  v;ere  studied  experimentally. 
The  static  and  dynamic  pressures  through  the  wake  were  ob- 
tained with  calibrated  Prandtl-type  pitot-static  probes, 
and  a  computer  program  was  written  to  extract  from  these 
data  the  velocity  profiles,  displacement  thickness,  momentum 
thickness,  and  the  shape  factor.   This  program  further  pro- 
vides the  turbulent  shear  coefficient,  eddy  viscosity,  and 
Prandtl  mixing  length  within  the  wake. 

A  wind  tunnel  was  designed  specifically  for  this  experi- 
ment and  its  construction  details  are  included.   The  tunnel 
test  section  is  unique  in  that  it  can  be  adjusted  to  generate 
any  specified  adverse  or  favorable  pressure  gradient. 

The  results  of  the  experiment  indicated  that  the  eddy 
viscosity  was  a  function  of  the  history  of  the  flow  rather 
than  of  the  local  mean  flov;  field,  as  is  classically  assumed. 
The  reaction  of  the  turbulent  wake  to  a  favorable  pressure 
gradient  after  encountering  an  adverse  pressure  gradient  was 
significantly  different  from  the  reaction  of  the  wake  to 
only  a  favorable  pressure  gradient. 
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NOMENCLATURE 

A    =  cross-sectional  area 

I  I   =  indicates  absolute  value 

A    =  turbulent  exchange  coefficient  (Eq.  3,  Introduction) 

a    =  constant  in  Prandtl's  eddy  viscosity  (Eq.  5,  Intro- 
duction) 

b    =  width  of  mixing  region  (Eq.  5,  Introduction) 

3    =  constant  in  test  section  area  distribution 
(Eq.  7,  Appendix  A) 

C,    =  volumetric  flow  rate  through  the  test  section 

Cp    =  negative  of  the  test  section  pressure  gradient 
(Appendix  A) 

1   2 
C    =  turbulent  shear  coefficient  =  x/p-pu  e 


d_ 

dx 


wake  generator  diameter 

indicates  differentiation  with  respect  to  x 


•sr—  =  indicates  partial  differentiation  with  respect  to  x 

dX 

6*  =  displacement  thickness 

e  =   eddy   viscosity   =   t/tt^ 

exp  =  indicates  exponentiation 

P  =  arbitrary  example  function 

H  =  shape  factor  =  6«/0 

i  =  indicates  i  direction  (i  =  1,2,3) 

»  =  large  distance  in  the  cross-streamwise  direction 

j  =  indicates  j  direction  (j  =  1,2,3) 

L  =  length  of  pressure  gradient 

Z  =  Prandtl's  mixing  length  (Eq.  4,  Introduction) 


In  =  logarithm  to  base  e 

y  =  laminar  coefficient  of  viscosity 

o  =  wake  centerline  conditions 

p  =  static  pressure 

1   2 

q  =  dynamic  pressure  =  p-pV 

p  =  fluid  density 

a.  .  =  stress  tensor  in  Navier-Stokes  equations 

a.  .  =  mean  flov/  stress  tensor  in  the  time  averaged 


ij 


Navier-Stokes  equations 


t  =  indicates  time 

T  =  indicates  the  characteristic  period  used  for 
time  averaging 

T..  =  Reynolds  shear  stress  tensor 

9  =  momentum  thickness 

u  =  mean  velocity  in  the  x-direction 

u'  =  fluctuating  velocity  in  the  x-direction 

u. ,u.  =  mean  velocity  in  i,j  directions  (i,j  =  1,2,3) 

ul,ul  =  fluctuating  velocity  in  i,j  directions  (i,j  =  1,2,3) 

V  =  mean  velocity  in  the  y-direction 

v'  =  fluctuating  velocity  in  the  y-directlon 

V  =  one  dimensional  mean  flow  velocity 
X  =  streamwise  coordinate 

X  =  j  coordinate 

y  =  cross-streamwise  coordinate 

y  =  lateral  position  at  v;hich  C   is  being  calculated 
(Eq.  5,  Appendix  B) 


Subscripts 

e  =  local  external  flov/  conditions 

l,j  =  directionality  parameter  (i,j  =  1,2,3) 

max  =  maximum  value  within  the  wake 

min  =  minimum  value  within  the  wake 

T  =  shear  related  quantity 

x,y  =  indicates  partial  differentiation  with  respect  to 
X  and  y 

1  =  test  section  inlet  conditions 

2  =  test  section  exit  conditions 

Superscripts 

=  indicates  a  time  averaged  quantity 
'     =  designates  a  fluctuating  quantity 


I.   INTRODUCTION 

When  tangential  stresses  are  applied  to  a  fluid  having 
internal  friction,  shearing  motions  are  set  up  in  line  with 
the  stresses  and  in  conformity  to  the  shape  of  the  bound- 
aries.  Within  this  flow  field,  various  kinds  of  secondary 
motions  become  possible.   Regular  ones,  when  they  occur, 
are  readily  accounted  for  in  terms  of  pressure  gradients 
resulting  from  the  curvature  of  the  main  flow,  usually 
imposed  by  the  shape  of  the  boundaries.   Irregular  ones, 
called  turbulence,  are  by  far  the  more  common,  and  their 
direct  cause  is  less  obvious.   Their  occurrence  does  not 
depend  on  the  shpae  of  the  boundaries,  but  like  all  secon- 
dary flov;s  they  must  depend  on  a  generating  mechanism  within 
the  flow  itself.   According  to  present  evidence  the  initial 
onset  of  turbulence  occurs  suddenly  by  a  breakdown  of  the 
laminar  flov;  in  localized  regions.   Turbulence  is  convected 
downstream  in  the  manner  of  any  other  fluid  property  and 
washed  av;ay  from  its  point  of  origin.   Repeated  breakdowns 
are  therefore  generally  required  to  maintain  a  continuous 
supply  of  turbulence,  and  instability  of  laminar  flow  is 
an  essential  part  of  this  process. 

In  treating  turbulence  mathematically  it  is  generally 
assumed  that  the  motion  can  be  separated  into  a  mean  flow, 
whose  components  are  u.  (Cartesian  tensor  notation),  and  a 
superimposed  fluctuating  motion  whose  components  are  u I ,  the 


mean  values  of  which  are  zero.   For  the  flow  of  an  incom- 
pressible Newtonian  fluid,  this  decomposition  can  be  applied 
to  the  dependent  variables  of  the  Navier-Stokes  equations 
and  time-averaged,  the  latter  defined  as 


t+T 
F(x,y)  E  ^  /   F(x,y,t)dt  (1) 

t 


where  T  is  a  period  of  time  large  compared  with  the  char- 
acteristic time  of  the  turbulent  oscillations.   The  resulting 
time-smoothed  Navier-Stokes  equations  are 


p5^=  ny  (°ij -p"i"j  -  p"l"j)         (2> 


These  equations  have  the  same  form  as  the  original 


Navier-Stokes  equations  if  a.,  is  replaced  by  a .  .  -  pu.'ul, 
Thus  the  equations  for  the  mean  flow  are  the  same  as  the 
original  equations  of  motion,  except  that  there  are  addi- 


tional virtual  stresses,  t  .  =  -pu.'ul,  which  represent  the 
mean  rate  of  transfer  of  momentum  across  a  surface  due  to 
the  velocity  fluctuations.   These  virtual  stresses  were 
first  introduced  by  Osborne  Reynolds  and  are  known  by  his 
name  [Ref.  4]. 

The  quantities  u!  and  ul  are  the  velocity  fluctuations 


in  the  i-  and  j-  directions  at  one  point  in  space  and  u'ul 
is  a  measure  of  the  correlation  between  the  fluctuations. 
Correlations  between  fluctuations  at  two  points  in  space 
are  also  measurable  and  give  the  size  and  orientation  of 
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the  turbulent  eddies.   In  trying  to  solve  analytically  for 
the  second-order  correlation  tensor,  the  solution  is  given 
in  terms  of  third-order  correlations.   Similarly,  the 
solution  for  the  third-order  correlation  is  given  in  terms 
of  the  fourth-order  correlation,  and  so  on.   Hence  the 
system  cannot  be  closed  and  it  is  impossible  to  solve 
exactly  for  the  second-order  correlations  [Ref.  1]. 

The  various  methods  for  the  calculation  of  the  turbulent 
flow  field  developed  to  date  are  based  on  empirical  hypoth- 
eses.  These  approaches  endeavor  to  establish  a  relationship 
between  the  Reynolds  stresses  produced  by  the  fluctuating 
velocities  and  the  mean  values  of  the  velocity  components 
and  their  spatial  derivatives.   By  this  means,  the  virtual 
stresses  are  given  an  analytical  form  which  closes  the 
system  of  equations  and  constitutes  the  starting  point  for 
the  calculation  of  the  mean  flow. 

The  first  hypothesis  relating  the  Reynolds  stress  to 
mean  flow  properties  was  advanced  by  Boussinesq  and  intro- 
duced a  turbulent  mixing  coefficient,  A  ,  such  that  the 
turbulent  shear  stress  could  be  written 


^  =  \  i  (3) 


This  coefficient  corresponds  to  the  viscosity,  y,  in  laminar 
flow  and  is  therefore  often  called  the  "virtual"  viscosity. 
Equation  (3)  has  the  great  disadvantage  that  the  virtual 
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viscosity,  A^,  is  not  a  property  of  the  fluid,  as  y  is,  but 
depends  on  the  structure  of  the  flow  field. 

In  1925  Prandtl  proposed  his  mixing  length  hypothesis 
in  which  the  turbulent  stresses  are  modelled  as  f ollov;s : 


^  =  p^'  if  if  w 


The  mixing  length,  £,  is  that  distance  in  the  transverse 
direction  which  must  be  covered  by  an  aggregation  of  fluid 
particles,  travelling  with  its  original  mean  velocity,  in 
order  to  make  the  difference  between  its  velocity  and  the 
velocity  in  the  new  lamina  equal  to  the  mean  transverse 
fluctuations  in  the  turbulent  flow.   The  mixing  length  is 
still  not  a  property  of  the  fluid,  but  it  is,  in  theory, 
a  purely  local  function  and  it  is  far  simpler  to  make 
plausible  assumptions  concerning  the  mixing  length,  Z^    than 
for  the  virtual  viscosity,  A  .   This  constitutes  the  essen- 
tial superiority  of  Prandtl 's  equation  over  that  due  to 
Boussinesq.   Even  so.  Equation  (^)  is  unsatisfactory  in  that 
the  apparent  kinematic  viscosity,  e  =  £  |-p-|  vanishes  at 
points  where  |-r^|  is  equal  to  zero,  that  is  at  points  of 
maximum  and  minimum  velocity,  which  is  contradictory  to 
experimental  evidence.   A  simpler  expression  for  e,  also 
proposed  by  Prandtl  which  overcomes  this  problem  is 


e  =  ab(u_^^  -  u^.  )  (5) 

max    mm 
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and  is  valid  only  for  free  turbulence,  where  a  is  an  experi- 
mentally determined,  dimensionless  constant  and  b  is  the  width 
of  the  mixing  zone.   Thus  e,  for  all  the  classical  models  of 
Reynolds  stresses,  remains  constant  over  the  v;hole  width  of 
the  wake  and  is  a  function  of  the  local  mean  flow. 

Several  other  similar  hypotheses  have  been  proposed  as 
analytical  expressions  for  the  Reynolds  stresses  which  will 
only  be  mentioned  here.   They  include:  Taylor's  vorticity 
transport  theory,  Von  Karman ' s  similarity  hypothesis,  and 
the  universal  velocity  distribution  law,  each  of  which  is 
linked  to  the  mixing  length  concept  due  to  Prandtl.   Each 
of  the  af orem.entioned  theories  is  valid  only  for  free  turbu- 
lent flow.   While  none  of  the  theories  adequately  models  the 
turbulent  process  in  general,  they  are  useful  tools  in  the 
integration  of  the  equations  of  motion  in  special  cases. 

In  this  investigation  it  was  desired  to  ascertain,  at 
least  qualitatively,  v:hether  the  virtual  viscosity  could  bf 
considered  to  be  a  function  of  the  local  mean  flow  propert  es 
in  the  case  of  two-dimensional  turbulent  wakes  in  constant 
favorable  and  adverse  pressure  gradients.   A  wind  tunnel 
was  designed  and  fabricated  for  this  study  which  incorporated 
a  variable-geometry  test  section  that  could  be  set  to  provide 
any  specified  pressure  gradient.   While  this  study  v;as  carried 
out  using  pitot-static  probes  as  data  sensors,  it  is  antici- 
pated that  the  flexibility  and  utility  of  the  wind  tunnel 
test  section  will  allow  a  very  detailed  analysis  of  the 
turbulent  wake  using  a  hot-wire  anemometer.   This  will  also 
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permit  the  assessment  of  more  modern  eddy  viscosity  models 
which  Involve  hypotheses  with  regard  to  the  dependence  of 
the  shear  on  the  turbulent  kinetic  energy. 
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II.   THE  WIND  TUNNEL 

A.   INLET 

The  inlet  bellmouth  shov/n  in  Fig.  1  was  constructed 
from  eight-inch  diameter  aluminum  pipe  halves  formed  into  a 
two-feet  square  inlet  section.   The  corner  seams,  which  were 
exposed  to  the  inlet  air  flow,  were  filled  with  epoxy  putty, 
hand  formed,  and  sanded  into  a  smooth  fillet.   This  was  done 
to  smoothly  transition  the  flow  around  these  corners  into 
parallel  flow. 

Immediately  downstream  of  the  inlet  bellmouth  was  the 
plenum  section  of  the  tunnel.   It  consisted  of  four  two-feet 
square  by  one  foot  long  steel  sections  bolted  together. 
Cork  gaskets  were  inserted  between  the  sections  to  prevent 
the  influx  of  air  at  the  joints  and  allow  for  smooth  air 
flow  over  these  junctions.   Finally,  two  copper  screens  of 
#16  mesh  between  the  first  two  and  last  two  sections  reduc  d 
the  ambient  turbulence  level  in  the  flov/  prior  to  its 
entering  the  nozzle. 

Initial  operation  of  the  tunnel  indicated  the  presence 
of  swirling  in  the  air  flow  through  the  test  section.   This 
swirling  was  eliminated  by  the  insertion  of  a  four-inch 
thick,  two-feet  square  section  of  quarter-inch  honeycomb 
just  prior  to  the  first  copper  screen. 
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B.   NOZZLES 

Two  different  nozzles  were  constructed  and  used  during 
the  course  of  this  investigation.   One  nozzle,  used  with  the 
favorable  pressure  gradient  test  section,  was  constructed 
from  half-inch  Plexiglas.   The  nozzle,  shown  in  Pig.  2,  is 
two  feet  long,  has  an  exit  area  of  l68  square  inches,  (24 
inches  high  by  7  inches  wide),  and  a  contraction  ratio  of 
3.^2.   The  top  and  bottom  sides  of  the  nozzle  are  parallel 
while  the  vertical  sides  contract  to  form  the  nozzle.   The 
wake  generator  was  attached  to  the  nozzle  in  its  exit  plane 
on  the  axial  centerline  of  the  vertical  walls.   Static 
pressure  ports  were  located  at  the  nozzle  entrance  and  exit 
planes  along  the  axial  centerline  of  the  bottom  wall.   These 
were  used  to  determine  the  nozzle  exit  (test  section  entrance) 
dynamic  pressure,  which  determines  the  magnitude  of  the  test 
section  pressure  gradient  (see  Appendix  A). 

The  second  nozzle  was  made  from  steel  and  was  used  witl 
an  adverse  pressure  gradient  test  section  configuration. 
This  nozzle,  shov/n  in  Fig.  1,  is  two  feet  long,  has  an  exit 
area  of  forty  nine  square  inches  (seven  inches  square)  and 
a  contraction  ratio  of  11.76.   Static  pressure  ports  were 
located  on  the  right  side  of  the  nozzle  at  the  inlet  and 
exit  for  the  purpose  of  determ.ining  the  exit  dynamic  pressure. 
The  wake  generator  vjas  not  mounted  in  the  nozzle  because  of 
the  difficulty  encountered  in  trying  to  center  and  align  it 
in  the  metal;  rather  in  this  case  it  was  mounted  just  at 
the  entrance  of  the  test  section. 
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C.  WAKE  GENERATORS 

The  v;ake  generators  v/ere  smooth  Plexiglas  cylinders 
spanning  the  width  of  the  tunnel.   A  one  inch  diameter 
cylinder  was  used  with  the  steel  nozzle  and  a  two  inch 
diameter  cylinder  v;as  used  with  the  Plexiglas  nozzle;  thus 
"equivalent"  turbulent  wakes  (see  Results)  were  generated 
for  each  test  section  configuration.   The  wake  velocity 
defect  and  turbulence  resulted  from  flow  separation  of  the 
turbulent  boundary  layers  on  the  cylinders. 

D.  TEST  SECTION 

As  shown  in  Fig.  3,  the  test  section  consisted  of  two 
parallel,  vertical  side  walls  separated  by  flexible  and 
adjustable  top  and  bottom  walls,  all  made  of  Plexiglas.   The 
flexible  v;alls  v;ere  made  from  one-eighth  inch  Plexiglas  and 
fastened  to  the  nozzle  with  piano  hinges.   Eight  adjustm.en- 
points,  each  with  two  threaded  rods,  allowed  these  walls  tc- 
be  positioned  to  any  desired  area  distribution  and  horizon 
tally  levelled  so  that  the  test  section  geometry  would  re- 
main rectangular  throughout.   Appendix  A  contains  the  devel- 
opment of  the  relationship  between  the  pressure  gradient 
and  area  distribution,  based  on  a  one-dimensional  analysis. 
Table  I  contains  the  dimensions  necessary  to  set  the  area 
profiles  and  form  the  two  pressure  distributions  used  in 
this  study. 

The  right  hand  vertical  wall  of  the  test  section  contained 
thirteen  slots  through  which  the  pitot-static  probes  were 


17 


inserted  into  the  flov/.   These  slots  were  nineteen  inches 
high,  three-eighths  inch  v;ide  and  spaced  five  inches  apart. 
A  teflon  sliding  seal  prevented  the  influx  of  air  through 
the  slot  in  which  the  probe  v/as  being  traversed.   A  filler 
plug  for  each  unused  slot  was  carefully  hand  milled  and 
polished  to  make  it  conform  perfectly  with  the  flow  side  of 
the  wall,  thus  eliminating  the  possibility  of  the  plugs 
protruding  into  the  flov/  and  inducing  separation.   A  thin 
aluminum  strip  vjedged  over  each  filler  plug  insured  a  tight 
fit. 

The  left  hand  vertical  v/all  of  the  test  section  contained 
three  rectangular  access  holes  equally  spaced  along  the 
axial  centerline.   These  enabled  the  operator  to  remove  and 
insert  the  wake  generators,  measure  the  probe  position 
accurately,  and  level  the  lower  adjustable  v;all.   Also,  along 
the  axial  centerline  of  this  wall  were  thirty-three  static 
pressure  ports  spaced  one-and-a-half  inches  apart.   A 
manometer  bank  connected  to  these  ports  gave  a  visual  indi 
cation  of  the  linearity  and  uniformity  of  the  pressure  dis- 
tribution.  The  manometer  also  served  to  indicate  any  regions 
of  separation,  as  did  several  rows  of  tufts  which  were 
attached  to  the  walls  of  the  adverse  pressure  gradient  test 
section. 

E.   DIFFUSER  AND  MOTOR 

The  test  section  was  followed  by  two  seven-inch  square, 
four  feet  long  steel  sections,  a  diffuser  section  and  a 


centrifugal  blower  driven  by  an  electric  motor  shown  in 
Figure  ^.   The  diffuser  was  attached  to  the  blower  by  a 
rubber  buffer  section  which  reduced  the  amount  of  bivration 
transmitted  from  the  motor  to  the  tunnel.   The  centrifugal 
blower,  operated  at  three  thousand  revolutions  per  minute, 
pulled  air  through  the  tunnel  and  exhausted  it  vertically 
through  a  muffler.   The  total  mass  flow  through  the  tunnel 
was  controlled  by  a  damper  mechanism  just  prior  to  the 
muffler. 

E.   DATA  SENSORS 

All  data  were  taken  with  two  calibrated  Prandtl-type 
pitot-static  probes  J  one  measuring  three-and-a-half  inches 
and  the  other  eighty-three  hundredths  inch,  stem  to  tip, 
thus  allowing  acquisition  of  data  betv;een  test  section  data 
stations.   Two  U-tube  manometers  v;ere  connected  to  the  probes 
so  that  direct  readout  of  local  static  and  dynamic  pressure 
was  obtained.   These  manometers,  shown  in  Fig.  5,  could  be 
inclined  at  any  angle  between  two  and  forty-five  degrees, 
providing  a  wide  range  of  sensitivity.   A  similar  manometer 
was  connected  to  the  nozzle  static  pressure  ports  so  that  the 
test  section  inlet  dynamic  pressure  could  be  continuously 
monitored. 

The  traversing  mechanism  used  to  position  the  pitot- 
static  probes  in  the  wake  was  hand-operated  and  completely 
detachable  from  the  tunnel.   The  probe  was  clamped  firmly  in 
its  movable  support  fixture  on  the  traverse  which  was  moved 
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vertically  by  a  screv;  drive  mechanism.   Probe  position  could 
be  determined  to  within  a  hundredth  of  an  inch  v/ith  the  aid 
of  a  ruled  indicator  on  the  traverse. 
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III.   EXPERIMENTAL  PROCEDURE 

A.   TEST  SECTION  PREPARATION 

The  test  section  v/as  configured  to  produce  both  a 
favorable-only  pressure  gradient  (Figure  6),  and  an  adverse- 
favorable  pressure  gradient  (Figure  3)  over  its  entire 
length.   For  ease  in  setting  the  flexible  walls  to  the  de- 
sired profile,  two  reference  lines  were  scribed  at  the  top 
and  bottom  of  the  right  hand  vertical  wall  of  the  test  sec- 
tion.  All  dimensions  given  in  Table  I  are  referenced  to 
these  lines  and  indicate  the  distance  in  inches  from  these 
reference  lines  to  the  flow  side  of  the  flexible  walls. 

The  following  procedure  was  used  in  setting  the  flexible 
walls.   Proceeding  from  the  test  section  inlet,  the  right 
side  of  each  adjustment  fixture  was  set  to  the  desired  posi- 
tion and  tightened.   After  all  eight  fixtures  had  been  set, 
each  was  re-checked  to  insure  that  its  setting  had  not  been 
disturbed.   Then,  with  the  aid  of  a  bubble  level,  the  left 
side  of  each  adjustment  fixture  was  moved  until  the  wall  was 
level.   Even  though  the  walls  could  only  be  set  precisely 
to  the  desired  dimensions  at  the  eight  adjustment  points, 
the  resulting  area  distribution  in  the  test  section  never 
varied  more  than  five  percent  from  the  desired  area.   The 
only  variance  occured  when  the  axial  area  gradient  was  large, 
such  as  at  the  beginning  of  the  favorable  pressure  gradient 
profile  or  in  the  transition  section  of  the  adverse-favorable 
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pressure  gradient  profile;  otherwise  the  wall  profile 
conformed  exactly  to  the  desired  area  distribution. 

B.   DATA  ACQUISITION 

All  data  were  recorded  from  the  water-filled  U-tube 
manometers  (Figure  5)  connected  to  the  pitot-static  tubes 
and  nozzles.   A  scale  was  placed  between  the  tubes  of  each 
manometer  providing  a  minimum  reading  of  a  hundredth  of  an 
inch  deflection  on  the  slant  height.   To  eliminate  any 
possible  parallax  and  to  assure  consistent  readings  from 
the  manometers,  a  portable  plastic  sight  unit  v;as  placed 
over  the  manometer  scale.   This  unit  restricted  one's  line 
of  sight  to  a  plane  perpendicular  to  the  manometer. 

At  each  position  of  the  probe  in  the  wake,  the  manometer 
readings  vjere  recorded  and  plotted  on  a  graph  so  that  it  was 
immediately  obvious  when  more  data  was  needed  in  a  particular 
portion  of  the  wake.   This  technique  was  particularly  advan- 
tageous in  determining  the  exact  center  of  the  wake  and  the 
point  where  the  freestream  was  encountered. 

The  high  turbulence  intensity  in  the  central  core  of  the 
wake  caused  large  fluctuations  in  the  manometer  water  columns, 
These  fluctuations  were  often  as  great  as  two-tenths  inch 
even  at  manometer  inclinations  of  thirty  degrees.   By   ob- 
serving the  fluctuations  for  several  minutes  at  each  probe 
position  in  the  wake,  it  was  possible  to  "average"  these 
fluctuations  by  eye  and  obtain  a  reading  that  was  accurate 
to  within  an  estimated  five-hundredths  of  an  inch. 
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The  inclination  (set  V7ith  an  inclincrr.eter)  cf  each  of 
the  three  U-tube  manor.eters  V7as  varied  according  to  tunnel 
operating  conditions  and  axial  location  of  the  pitot-static 
probe.   For  a  given  axial  data  station  and  local  flow  veloc- 
ity, an  inclination  was  chosen  that  would  yield  approximately 
a  ten  to  eleven  inch  difference  in  vrater  colurji  slant  height 
when  the  probe  was  positioned  in  the  freestream.   This 
assured  adequate  sensitivity  and  readability  of  the  vjater 
columns  at  the  lov:er  velocities  cf  the  v.'ake  centerline. 

Traversing  of  the  v/ake  by  the  pitot-static  probe  was 
alv;ays  done  fror.  the  upper  free  stream  edge  to  the  lower 
free  stream  edge.   The  free  stream  was  located  by  moving 
the  probe  away  from  the  upper  v:all  until  the  maxir.uri  slant 
height  difference  on  the  dynaFiic  pressure  r.anor^eter  remained 
unchanged  for  a  fev:  tenths  of  an  inch  belov;  the  point  at 
which  it  first  occurred.   This  point  v;as  taken  as  the  free 
stream.   The  lower  free  stream  v;as  found  in  the  sar.e  way. 
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IV.   DISCUSSION  OF  RESULTS 

A.  TEST  SECTION 

A  qualitative  indication  of  the  linearity  of  the  test 
section  pressure  gradient  in  both  the  favorable  and  the 
adverse-favorable  configurations  can  be  seen  in  Figures  7, 
8.   The  limitations  of  the  test  section  in  generating  an 
exactly  linear  pressure  distribution  result  from  the  lack 
of  completely  flexible  upper  and  lower  test  section  wall 
and  from  the  limited  number  of  wall  adjustment  points  from 
which  to  set  the  area  profile.   Considering  these  physical 
limitations  and  the  fact  that  the  desired  area  distribution 
in  the  test  section  v/as  derived  from  one-dimensional  analysis, 
the  resulting  actual  pressure  distribution  in  the  test 
section  is  quite  good.   A  comprehensive  program  to  calibrate 
the  test  section  was  not  undertaken  since  only  a  qualitative 
idea  of  its  performance  was  necessary  for  this  investigation. 

B.  EXPERIMENT 

The  experiment  was  designed  to  compare  the  decay  of  two 
similar  wakes  which  have  had  different  flow  histories.   This 
was  done,  within  the  physical  constraints  of  the  test  section, 
by  using  a  two-inch  diameter  wake  generator  in  a  favorable 
pressure  gradient  channel  sixty  eight  inches  long  and  a  one- 
inch  diameter  wake  generator  in  an  adverse-favorable  pressure 
gradient  channel  with  the  favorable  pressure  gradient 
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occupying  the  downstream  thirty  four  inches  (Figures  3  and 
8).   The  two  wakes  are  compared  within  the  favorable 
pressure  gradient  regions  downstream  of  the  point  where 
both  wakes  have  the  same  velocity  defect.   The  approximate 
solution  of  the  integral  momentum  equation  in  Appendix  C 
illustrates  the  possibility  that  an  adverse  pressure  gradient 
can  sustain  a  v;ake  deficit,  such  as  the  one  studied  here. 
In  addition,  the  v/akes  are  compared  at  the  same  relative  non- 
dimensional  distance  into  the  favorable  pressure  gradient. 
This  was  possible  because  the  favorable  pressure  gradient 
lengths  were  in  the  same  ratio  as  the  initial  wake  widths. 
It  is  apparent  from  the  experimental  data  in  Figure  9  that 
the  velocity  defect  decays  rapidly  v/hen  the  flovj  enters  the 
favorable  pressure  gradient  directly  from  the  nozzle.   How- 
ever, the  velocity  defect  persists  [Figure  10]  and  then 
decays  at  a  much  slov;er  rate  [Figure  11]  when  the  flow  has 
previously  encountered  an  adverse  pressure  gradient.   Sine. 
the  rate  at  which  the  velocity  defect  diminishes  is  direct  .y 
related  to  the  eddy  viscosity,  it  can  be  at  least  qualita- 
tively concluded  from  Figures  12,13  that  the  eddy  viscosity 
of  a  turbulent  wake  is  a  function  of  the  flow  history  and 
its  effect  on  the  turbulent  field.   Classically,,  it  is 
assumed  that  the  eddy  viscosity  is  purely  a  function  of  the 
local  mean  flow  properties,  we  must  conclude  that  care  must 
be  exercised  when  applying  classical  eddy  viscosity  models 
to  pressure  gradient  situations.   In  order  to  test  this 


25 


conclusion,  a  numerical  analysis  of  the  experimental  data 
was  conducted. 

C.   COMPUTER  PROGRAM 

The  data  analysis  computer  program  uses  a  forward 
difference  scheme  to  solve  the  equations  of  motion  for  the 
turbulent  shear  coefficients  (CTAU)  and  the  eddy  viscosi- 
ties (EDDY)  as  described  in  Appendices  B  and  D.   The  exper- 
imental data  obtained  from  the  wind  tunnel  were  used  to 
test  the  program  and  assure  correct  processing  with  no 
coding  errors.   The  numbers  generated  for  CTAU  and  EDDY  were 
unreasonably  large  and  therefore  the  results  are  inconclusive 
It  is  felt  that  the  increments  in  the  data  in  the  axial 
direction  used  in  finite  differencing  the  data  were  too 
large  for  an  accurate  approximation  of  the  partial  deriva- 
tives which  appear  in  the  equations  (see  Appendix  B).   A 
better  method  would  be  to  smooth  the  data  and  use  a  central 
difference  numerical  scheme  to  approximate  the  partial 
derivatives.   The  data  taken  in  the  experimental  phase  of 
this  investigation  were  not  adaptable  to  the  central 
difference  scheme.   Given  experimental  data  with  sufficiently 
small  increments  in  the  axial  direction,  the  program  will 
provide  the  correct  values  for  the  parameters  listed  in 
Appendix  D. 
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V.   CONCLUSIONS 

The  experiment  described  in  this  report  was  designed  to 
test  the  validity  of  classical  eddy  viscosity  models  in 
free  turbulent  mixing  situations  with  axial  pressure  gra- 
dients.  The  results  indicate  that  those  classical  models 
which  are  explicitly  independent  of  the  local  turbulence 
structure  will  fail  to  predict  the  experimentally  observed 
phenomena  which  are  a  strong  function  of  the  turbulent  field. 
In  these  pressure  gradient  situations,  the  local  turbulent 
kinetic  energy  is  not  a  unique  function  of  the  local  non- 
turbulent  mean  flov;  and  the  influence  of  the  turbulence 
on  the  form  of  the  Reynolds  stresses  must  be  an  explicit 
one  in  any  eddy  viscosity  model. 
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Figure  1 
View  of  tunnel  showing  nozzles  used  with  (1)  fav- 
orable pressure  gradient  and  (2)  adverse-favorable 
pressure  gradient  test  section  configurations. 
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Figure  2 
View  of  Plexiglass  nozzle  and  test  section  showing 
(1)  2-inch  wake  generator  installed  in  the  nozzle 
and  (2)  the  pitot  tube  traverse  mechanism. 
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Figure  3 
View  of  Plexiglass  test  section  showing  (1)  the 
adjustable  walls  in  the  adverse-favorable  pressure 
gradient  configuration,  (2)  the  adjusting  rods  for 
these  walls  and  (3)  the  slots  through  which  the  pitot 
tubes  are  inserted. 
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Figure  k 
Overall  view  of  the  wind  tunnel  showing  the  diffuser, 
centrifugal  blower  and  its  muffler  at  extreme  right 
in  picture. 
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Figure  5 
View  of  Inclined  U-tube  manometers  used  to  obtain 
pressure  readings  in  the  nozzle  and  test  section. 
The  sight  unit  (1)  used  to  eliminate  parallax  in 
reading  the  manometers  is  shown  on  the  left. 
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Figure  6 
Test  section  configured  for  the  favorable  pressure 
gradient.   The  milled  slot  through  which  the  pitot 
tube  is  inserted  can  be  seen  at  station  ^. 
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Figure  7 
View  of  manometer  bank  connected  to  the  static 
pressure  ports  along  the  test  section  centerline 
Shown  is  the  linear  pressure  distribution  in  the 
test  section  resulting  from  the  configuration  in 
Pig.  6.   The  horizontal  water  columns  are  not 
connected  to  the  tunnel. 
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Figure  8 
View  of  manometer  bank  connected  to  the  static 
pressure  ports  along  the  test  section  centerline. 
Shown  is  the  pressure  distribution  in  the  test 
section  resulting  from  the  configuration  in  Fig.  3 
The  horizontal  water  columns  are  not  connected  to 
the  tunnel. 
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APPENDIX  A 
DESIGN  OF  THE  TEST  SECTION  AREA  DISTRIBUTION 

The  local  cross-sectional  area  of  the  test  section  and 
the  static  pressure  of  the  fluid  flov/ing  through  the  section 
are  related  through  the  equation  of  continuity  and  Bernoulli's 
equation.   For  the  steady,  isentropic  flow  of  an  incompres- 
sible fluid,  these  equations  can  be  written: 

VA  =  constant  (1) 


1   2 

p  +  p-pV  =  constant  (2) 


In  the  log-derivative  plane  (1)  and  (2)  become 


dV  _    dA 
V    "A 


(3) 


dp  =  -pV^  ^  (^) 


Substitution  of  (3)  into  (^)  yields 


,2  dA 


dp  =  pV^  ^  (5) 


From  (1)  V  =  CVA,  v;here  C^  is  the  volumetric  flow  rate  at 
the  inlet  to  the  test  section.   Thus,  (5)  can  be  written 


il2 


c  2 


After  specifying  a  particular  pressure  gradient,  (6)  can 
be  solved  to  determine  the  required  area  distribution. 

For  a  constant  pressure  gradient,  ^2.  =  _c   and  inte- 
gration of  (6)  for  this  simple  case  is  straightforward. 
The  result  can  be  written 


^1       1  +  3x 


where  3  =  Cp/q  ,  A,  is  the  test  section  inlet  area  and  q 
is  the  dynamic  pressure  at  the  test  section  inlet.   For  a 
specified  inlet  and  exit  area 

B  =  ~ p-^—  (8) 

A2   L 

where  L  is  the  length  of  the  pressure  gradient,  and  A^  is 
the  channel  area  at  x  =  L. 

The  area  profiles  for  the  favorable  and  adverse-favorable 
pressure  gradients  used  during  this  investigation  are  shown 
in  Figs.  3  and  6,  respectively.   Table  I  contains  the  specific 
dimensions  used  to  set  the  flexible  walls  of  the  test  section 
to  the  desired  area  profile.   Shown  for  comparison  in  this 
table  are  measurements  of  the  upper  and  lower  flexible  walls 
as  they  were  set  prior  to  taking  data.   The  differences  in 


^3 


dimensions,  where  they  exist,  between  the  theoretical  and 
actual  wall  settings  show  the  degree  of  conformity  of  the 
flexible  walls  to  the  desired  profile. 
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APPENDIX  B 

DEVELOPMENT  OP  AN  EQUATION  FOR  CALCULATION 
OP  THE  LOCAL  TURBULENT  SHEAR  COEFFICIENT 


The  momentum  equation  governing  the  mean  flow  properties 
of  a  steady,  tv70-dimensional,  incompressible  turbulent 
shear  flow  can  be  v;ritten,  using  Reynolds'  averaging 
technique  [Ref.  3]: 


^  9x   ^  8y     p  dx    p  8y  ^^^ 


where  t  =   -pu'v'   and  pu'   has  been  neglected  with  respect 
to  p.   The  same  arguments  v;hich  Prandtl  used  to  derive 
equation  (1)  from  the  Navier-Stokes  equations,  in  the  case 
of  a  boundary  layer,  are  directly  applicable  to  a  two- 
dimensional  wake  in  approximately  parallel  flow  through  a 
rectangular  duct.-  Introducing  the  integrated  form  of  the 
continuity  equation 


y   =  "  /(|^)  dy  (2) 

0    ^^ 


into  equation  (1)  transforms  it  to 


-  au  _  9u   /  /3U)  d^  =  _  1  dp   +  1  9t      ,3) 
"^  87   8y  Q-'   ^3x^  "^y     p  dx  ^  p  9y     ^^^ 


By  non-dimensionalizing  equation  (3)  with  respect  to  the 
local  freestream  velocity,  u^,  and  utilizing  the  product  rule 
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of  differential  calculus,  equation  (3)  becomes 


u  u 

e      —     V  —        —    ( 

e  e  e   ^    0  e  e  *^   e   0    u 


U    /  U  X     ,   ,  U  n2    ^X      /Us       f^/  U  V   .        /  U  N   /   ^Xn      ^z  - 


e 

F    +  -TT^  (1) 


2 
where  C   =  2t/pu   and  subscripts  x  and  y  indicate  partial 

differentiation  with  respect  to  x  and  y.   Equation  (k)    is 
now  integrated  from  y  =  0  to  y  =  y  ,  an  arbitrary  point  in 
the  wake,  and  relax  the  condition  that  equation  (M  be 
satisfied  at  every  point  in  the  flow,  allowing  it  to  be 
satisfied  on  the  average  between  y  =  0  and  y  =  y  .   After 
the  necessary  integration  by  parts  and  considerable  manipu- 
lation, the  final  equation  for  the  turbulent  shear  coeffi- 
cient, C  ,  can  be  written 

T 


^.^^^  =  s\,^,   .  5ll3C(^).  .  (^)(%]dy  .  /^-l^F  dy 


2     0     ^e    ^e    ^e  ^   ^e  ^e      0   pu/  ^ 

(5) 

Equation  (5)  was  solved  numerically  with  the  computer 
program  listed  in  Appendix  D. 

The  validity  of  equation  (5)  can  be  checked  by  letting 


^T  - 
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Thus  C^(y^)  =  0  and  equation  (5)  reduces  to 


,Q    -,  du 
e 


where 


e  E   /  -IL  (1  _  JL)  dy 
0  %      ^e 


6*  E   /  (1  -  Ji.)  dy 
0       ^e 


Equation  (6)  is  the  well  known  momentum  integral  equation 
in  boundary  layer  theory. 
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APPENDIX  C 

SOLUTION  TO  THE  MOMENTUM 
INTEGRAL  EQUATION 


The  momentum  integral  equation  obtained  in  Appendix  B 
can  be  integrated  in  the  special  case  of  a  constant  pressure 
gradient.   Prom  Eq.  (6)  in  Appendix  B 


e 


where  6,  6*  are  the  momentum  and  displacement  thickness 

defined  in  Appendix  B. 

du 

An  expression  for  -r-^  can  be  obtained  from  the  inviscid 
^  dx 

Bernoulli  equation 


dx      pu^   dx 


e 


Substitution  of  Eq.  (2)  into  Eq.  (1)  yields  a  differential 
equation  in  terms  of  6,  ^,  and  H  =  (S*/e  (the  shape  factor) 


|_(ln  e)  -  -1^  (|£)  (H  +  2)  =  0       (3) 
e 


The  solution  of  Eq.  (3)  is 


/-  =  exp[  /  -^  (||)  (H  +  2)dx]       (4) 
1       0   pu^ 
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Prom  Appendix  A,  ^  =  -6q   and  the  integrand  of  Eq.  (^) 
can  be  written 


3q. 
(jT-^    (H  +  2)  (5) 

e 


From  the  results  of  Appendix  A,  the  ratio  is  a  known  function 
of  X 


Qg    1  +  Bx 


(6) 


After  substituting  Eqs.  (5),  (6)  into  Eq.  (^)  and  integrating, 
the  final  solution  becomes 


T-h^'-P^/      2(i'h-%x)   ''^^   (7) 


Thus,  experimentally  derived  values  for  the  shape  factor 
can  be  substituted  into  Eq.  (7),  and  the  resulting  value  of 
0  can  be  compared  v/ith  the  experimentally  obtained  values. 
For  the  restrictive  case,  H  =  constant,  we  can  extend 
this  solution  to 


l+H/2  ^^^ 


"1     (1  +  6x 

Further,  for  H  -  1.0  (a  close  approximation  under  the 
conditions  of  this  investigation),  Eq.  (7)  reduces  to 


^9 


(1  +  Bx)3/^ 


(9) 


This  result  compares  favorable  with  experimental  data  and 
serves  as  a  guide  for  calculating  the  growth  or  decay  of 
wakes  in  pressure  gradients. 
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APPENDIX  D 
DATA  REDUCTION  PROGRAM 

A.  GENERAL 

The  following  FORTRAN  IV  computer  program  uses  the 
experimental  values  of  dynamic  and  static  pressure  through 
the  wake  in  inches  of  v/ater  and  calculates  the  wake  momentum 
thickness,  displacement  thickness,  turbulent  shear  coeffi- 
cient, eddy  viscosity,  and  Prandtl  mixing  length.   The 
results  are  printed  out  as  well  as  plotted  on  the  calcomp 
plotter. 

B.  NOMENCLATURE 

The  following  parameters  and  arrays  are  used  in  the 

program,  with  their  respective  dimensions  shown  in 

parentheses : 

ALV    -  array  used  to  store  interpolated  values  of 
dynamic  pressures  (feet) 

AMIXL   -  local  Prandtl  mixing  length  in  the  wake  (feet) 

CTAU   -  local  turbulent  shear  coefficient  for  each 
point  in  the  wake 

DELRAT  -  ratio  of  succeeding  DELSTR's  to  that  at  the 
first  data  station 

DELSTR  -  wake  displacement  thickness  (feet) 

DELX   -  axial  distance  between  succeeding  probe  tip 
positions  (feet) 

DIAMTR  -  diameter  of  wake  generator  (inches) 

EDDY   -  eddy  viscosity 


51 


EPSIL  -  residual  in  momentum  integral  equation  when 
experimental  values  for  THETA  and  SHAPE  are 
substituted  into  it 

LABEL   -  label  for  subroutine  DRAW 

MN  -  number  of  data  points  laterally  through  the 
wake 

MX     -  number  of  axial  data  positions 

MY     -  number  of  points  (always  odd)  into  which  the 
wake  is  divided 

P      -  array  of  local  static  pressures  for  each  point 
in  the  wake  (psf) 

PATM   -  atmospheric  pressure  (inches  of  mercury) 

Q      -  freestream  dynamic  pressure  (psf) 

SHAPE   -  wake  shape  factor  (DELSTR/THETA) 

THETA   -  wake  momentum  thickness  (feet) 

THETAP  -  inclination  of  static  pressure  manometer  (degrees) 

THETAV  -  inclination  of  dynamic  pressure  manometer 
(degrees ) 

THRAT   -  ratio  of  succeeding  THETA' s  to  that  at  the 
first  data  station 

TITLEl  -  title  for  subroutine  DRAW 

TITLE4  -  title  for  subroutine  DRAW 

U      -  array  of  velocities  for  each  point  in  the 
wake  (feet  per  second) 

URATIO  -  ratio  of  local  wake  velocity  to  local  free- 
stream  velocity 

XIN    -  array  of  probe  traverse  positions  through  the 
wake  (inches) 

XOUT   -  interpolated  probe  traverse  positions  for  each 
of  the  MY  points  in  the  wake  (inches) 

XPOS  -  probe  tip  position  measured  from  downstream 
side  of  the  wake  generator  (inches) 
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YOUT   -  interpolated  manometer  readings  (for  YQ  and 
YPS)  for  each  of  the  MY  points  through  the 
wake  (inches) 

YPS    -  static  pressure  manometer  slant  distance 
height  difference  (inches) 

YQ     -  dynamic  pressure  manometer  slant  distance 
height  difference  (inches) 


C.   SUBROUTINES 

INTEG  -  uses  Simpson's  rule  to  numerically  integrate 
the  function  F.   It  must  be  supplied  with  K 
points  (alv/ays  odd)  equally  spaced,  DELX 
units  apart. 

PIF2    -  second-order  interpolation  function.   X  is  the 
point  for  which  interpolation  is  desired. 
XLIST  and  FLIST  contain  the  N  number  of  abscissas 
and  ordinateSj  respectively,  of  the  table  to 
be  interpolated. 

DECIDE  -  insures  that  the  wake  is  divided  into  an  odd 
number  of  points  spaced  0.1  inches  apart.   If 
the  wake  contains  an  even  num.ber  of  experimental 
data  points  (found  by  subtracting  the  lower 
traverse  position  from  the  upper  position  and 
multiplying  by  ten),  DECIDE  adds  a  point  onto 
the  upper  edge  of  the  wake  and  sets  its  value 
equal  to  the  freestream  value.   The  total 
number  of  points  are  printed  out. 


FILL 


DRAW 


insures  that  th 
data  positions 
the  one  before 
precluded  in  th 
EPSIL  and  CTAU. 
has  fewer  point 
uppermost  wake 
into  additional 
number  of  point 


e  wake  at  succeeding  downstream 
contains  as  many  data  points  as 
it.   Thus  division  by  zero  is 
e  computation  scheme  used  for 

If  the  succeeding  data  station 
s  in  its  wake,  FILL  copies  the 
point  (the  freestream  point) 

points  as  necessary  to  equal  the 
s  at  the  current  position. 


SSP3  subroutine  for  the  IBM  360/6?  at  the  Naval 
Postgraduate  School  which  uses  the  CALCOMP 
plotter  to  graph  the  program  output. 
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D.   USAGE 

The  program  villi   accept  up  to  twelve  data  positions  in 
the  streamwise  direction.   If  more  than  that  number  are 
required,  the  dimensions  of  the  appropriate  arrays  (ALP, 
ALV,U,P,A,B,C,D,EDDY,AKIXL,XOUT,URATIO,CTAU,DELSTR,THETA, 
SHAPE , EPSIL , Q , DELRAT , THRAT , THETAV , THETAP , P ATM , XPOS , DELY , 
NUM)  will  have  to  be  modified  accordingly.   The  program 
is  set  up  for  exactly  three  complete  runs  (a  run  is  defined 
as  the  processing  of  data  taken  for  a  particular  test 
section  configuration  and  inlet  dynamic  pressure).   If  more 
or  fev/er  runs  are  necessary,  the  limits  for  the  KK  DO-LOOP 
have  to  be  modified  accordingly.   There  is  no  restriction 
on  the  number  of  data  points  that  can  be  obtained  at  each 
axial  data  position.   Subscript  I  alv/ays  refers  to  each 
axial  position  for  v/hich  data  was  taken,  and  subscript  J 
always  refers  to  the  points  through  the  wake,  at  a  given 
axial  data  position. 

The  following  input  data  is  supplied  by  the  user  in  the 
order  presented.   Input  format  is  given  in  parentheses. 

Two  cards:   TITLEl  (6a8)  -  first  graph  title 

Two  cards:   TITLE4  (6a8)  -  second  graph  title 

One  card:    LABEL  (12A4)  -  curve  label.   One  number  in  each 
four  colujTins  to  designate  the  data  position, 
i.e.   12   3^,  etc. 

For  each  run  the  following  data  is  needed: 

One  card:    run  identification  (up  to  Hj   characters) 
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One  card:    number  of  data  positions  (12);  wake  generator 
diameter,  DIAMTR  'Inches  ,  (F10.2) 

One  or  more  cards:  THETAV,  degrees  (8F10.1) 

One  or  more  cards:  THETAP,  degrees  (8F10.1) 

One  or  more  cards:  PATM,  Inches  (8F10.2) 

One  or  more  cards:  XPOS,  Inches  (8F10.2) 

One  card:    data  position  Identification  (up  to  ^7  characters) 

One  or  more  cards:  XIN,  Inches  (8F10.3) 

One  or  more  cards:  YQ,  Inches  (8F10.3) 

One  or  more  cares:  YPS,  Inches  (8F10.3) 

E.   PROGRAM  LISTING 

A  listing  of  the  program  v/hlch  has  been  described  here 
is  Included  on  the  follov;ing  six  pages. 
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IMPLICIT    REAL*4     (A-HtO-Z) 

REALMS       TITLE1(12) 

REAL*8       TITLE4(12) 

REAL*8        LABELl/'  DELS'/ 

REALMS       LABEL2/'  THET'/ 

REAL'i'8   LABEL3/'     SHAP«/ 

REAL*4       LABELdZ) 

DIMENSION    ALP{ 12,125) tALV( 12,125) ,U( 12,125) ,P{ 12, 125) 

DIMENSION    A(12,125),B{12,i25),C(12,I25),D{12,12  5) 

DIMENSION    EDDY (12,125)  , AM  IX L ( 12 , 125 ) , XOUT ( 12, 125) 

DIMENSION    URATI0(12,125) ,CTAU( 12,125),YQ(30) ,YPS(30) 

DIMENSION  UR(125), UPC  125), UC( 125), UD{125),F(1 25) 

DIMENSION  Y0UT{125)  ,XI  NOG) 

DIMENSION  AXIS( 125) 

DIMENSION  DELSTR(12),THETA(12),SHAPE{12),EPSIL(12) 

DIMENSION  0( 12) ,DELRAT(12) ,THRAT(12) ,THETAV{12) 

DIMENSION  THETAP (12),PATM{ 12) ,XPOS{ 12) ,DELY(12) 

DIMENSION  NUM(12) 

COMMON/A/  H 

EQUIVALENCE  { F ( 1  )  , YOUT ( 1 ) , UR { 1 ) , UP ( 1 ) , UC ( 1 ) , UD  (  1)  ) 

EQUIVALENCE   ( A( 1 ,1) , EDDY (1,1) ) 

EQUIVALENCE   ( B (  1 , 1 ) , A M 1  XL ( 1 , 1 )  ) 

EQUIVALENCE  (C  ( 1 ,  1  )  ,  AX  IS ( 1 ) ) 

DATA  GAMA/62. 32/, RHO/0. 002378/ 

H=0.1E0 

Y=H/12.0 

50  FORMAT (• 
1'  ) 

51  FORMAT( I2,F10.2) 

52  FORMAT (BFlO.l ) 

53  FORMATC 8F10.2) 

54  FORM AT (« 
1'  ) 

55  F0RMAT(I2) 

56  FORMAT(8F10.3) 

57  FORMAT(///////  '      •,T20,6A8) 

58  FCRMAT( /////// •     ',T20,12A4) 

70  FORMAT(//////////'      •  , T33 , • ECHO-CHEC K    DATA • /// ) 

71  FORMAT(//'     •  ,T20,'THETAV'  ,  ('      •  .  T33 , 8F 1  J . 2/ / ) ) 

72  FORMAT(//'     ', T20 ,' THETAP ',( •      • , T33 , 8F 1 0. 2// ) ) 

73  FORMAT(//'     • ,T 20 , • P ATM • , ( •      ' , T33, 8F 10. 2/ / ) ) 

74  FORMAT(//«     ' , T2 0 , • X POS • , ( •     • ,T33, 8F1 0 . 2/ / ) ) 

75  FORMAT(///»     • , T53 ,  ' MN= •  ,  I  3///) 

76  FORMATC//'     •  ,T 1 0 , • X  IN • ,  ( '     « , T2 1 , 8F1 0. 3// ) ) 

77  FORMATC//'     ',T10,'Y0    ',('     ' ,T2 1 , 8F1 0.3/ / ) ) 

78  FORMATC//'     ' , T 10 , ' YPS ' , C '      ' , T2 1 , 8F 1 0. 3/ / ) ) 

79  FORMATC  •  1«  ,T50, '*>.-*RESULTS    OF    DATA    RE  DUCT  ION=!^*v  •  ) 

80  FORMATC //////•     ',T50, 'AXIAL    STATI ON '  , I  2 , ' OUT    OF '  , I  2 , 
1  'STATIONS'///// ) 

81  FORMATC     '  ,T18  ,  •  DEL  ST  AR='  ,  G12  .  5  ,  5X  ,  '  T  HET  A=  '  ,  G  12.  5,  5X, 
1 'SHAPE=' ,G12.5, 5X, ' EPSILON=' ,G12.5) 

82  FORMATC//////'      '  , T8 ,  ' PO INT ' , T 30 ,  ' U '  ,  T48  ,  ' U /UE ' , T70, • P  ' 
1  ,T88,'CTAU'  ) 

83  FORMATC  '0'  ,T9 , I  3 ,T 24,  G 12. 5 , T44, G 12.  5 , T64 , Gl 2 . 5 , T84 , 
1G12.5) 

84  FORMATC  •  1'  ,T44,  ":^**SUMMARY    OF     INTEGRATED    WAKE    PROPER', 
I'TIES*^?^'  ) 

85  FORMATC/////'     •  , T3 ,  ' AX  I AL ' , T75,  ' THET A • / '      ', T2 ,' STATI ON 
1'  ,T16,  'X/D'  ,T2  9, 'DELSTAR'  ,T44, • DEL RAT  10'  ,T60,' THET A'  , 
2T75, 'RATIO' ,T9  0, 'SHAPE' ,T105,'EPSILON') 

86  FORMATC '0',T4,I3,T11,G12.5,T26,G12.5,T41,G12.5,T56, 
1G12.5,T71,G12. 5,T86,G12.5,T102,G12.5) 

87  FORMAT  C  '  1*  ,T49,  '  ^-*'^REYNOLD"S    MIXING    LENGTH    DATA=i'**') 

88  FORMATC//////'      '  ,T6 1 , ' ST  AT  ION*  ,  13///  •     '  ,  T44,  ' PO I  NT • , 
1T59,'EDDY    V  I SC 0  SI  TV'  , T79 , ' Ml  XI NG    LENGTH') 

89  FORMATC     '  ,  T45  ,  13  ,  T60  ,  G  12.  5,  T  80  ,  G  12.  5) 

90  F0RMATC6A8) 

91  FORMATC//'  ' , ' LAST= • , I  3) 

92  FORMATC 12 A4) 


56 


READI5,90)  TITLEl 
READ(5,90)  TITLE4 
READ(5,92)        LABEL 

ZERO  ALL  ARRAYS  BEFORE  COMPUTATION  BEGINS 
FC J)=Z 


DO  1000 
Z  =  0.0 

DO  i  J  = 
DO  1  1  = 
XOUT{ I , 
ALP( I, J 
ALV(  I,  J 
U(I . J) = 
URATIO( 
P(I, J)= 
A(I  .  J)  = 
B{  I  ,J)  = 
C(I,  J)  = 
D(  I, J)  = 
CTAUd  , 
DELSTR( 
THETAf  1 
SHAPE(  I 
EPSIL(  I 
0( I)=Z 
DELRAT( 
THRAT(  I 
THETAVt 
THETAP{ 
PATM(  I  ) 
XPOS{  I  ) 
DELYd  ) 
CCNTINU 
DO  2  1  = 
XIN£  I}  = 
YQ(I  )  =  Z 
YPSd  )  = 
CCNTINU 
READ(5, 
READ(5, 
READ( 5, 
READ(5, 
READ(5, 
READ( 5, 
WRITE (6 
WRITE(6 
WRITE{6 
WRITE  (6 
WRITE (6 
WRITE(6 
WRITE{6 
WRITE(6 
WRITE(6 


KK=1,3 

ltl25 

1,12 

J)=Z 

)=Z 

)=Z 

Z 

I  ,J)=Z 

z 
z 
z 
z 

z 

J)=Z 
I  J=Z 
}  =  Z 
J=Z 
)  =  Z 

I  )=Z 
)=Z 
I  )  =  Z 
I  )=Z 

=  z 
=  z 
=z 

E 

1  ,30 

z 


50J 
51) 
52) 
52) 
53) 
53) 
,50 
,70 
,71 
,72 
,73 
,74 
,57 
,  57 
,58 


MX, DIAMTR 
(THETAV( I) ,I=1»MX) 
(THETAP{  I  )  ,  1  =  1, MX) 
(PATM(  I  )  ,1=1, MX) 
(XPOS(  I  )  ,1=1, MX) 


(THETAVd  )  ,1  =  1, MX) 
(THETAPt I  )  ,  1  =  1, MX) 
(  PATMd)  ,1=1,  MX) 
{ XPCSl  I  ),  I  =  1,MX) 
(TITLEK  I  )  ,1=1,6) 
{  TITLE4( I)  ,1=1,6) 
(LABEL! I ), 1=1, 12) 


BEGIN  COMPUTATION  FOR  EACH  STATION 

MX1=MX-1 

DO    6    1=1, MX 

THETAV{  I  )  =  THETAV( I  )/57  .296 

THE  TAP  (  I  )=ThETAPd  )/5  7.296 

PATM(  I  )=PATf',(I  )^70.73 

XPOS( I)=XPOS (I )/12.0 

READ    IN    MANOMETER    HEIGHTS 


READ(5,54) 
REAC(5, 55) 
READ(5,56) 
READ(5, 56) 
READ(5,56) 


MM 

(X  IN(J ) ,J  =  1,MN) 

( Y0( J)  ,J  =  1,MN) 

(YPS(J )  ,J  =  1,MN) 
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C      ECHO-CHECK  INPUT 
C 

WRITE(6,54) 

WRITE(6t75)   MN 

WRITE(6*76)   (XIN( J) ,J=1,MN) 

WRITE{6.77)   (YOCJ) ,J=ltMN) 

WRITE(6,78)   (YPS( J) ,J=1»MN) 

CALL     DECIDE(XIN,MNtMY) 

NUMd)  =MY 
C 

C  INTERPOLATE    FDR    MISSING    DATA    POINTS    INSIDE    WAKE 

C 

XOUTd  t  1J=XIN(  1) 

DO    21    J =2, MY 

XOUTd  »  J)=XOUT(  la  )+FLOAT(  J-1)*H 
21    CONTINUE 

DO  3  J=1,MY 

X=XOUT(  I  ,J) 

YOUT  (  J  )  =  P I  F  2  (  X  ♦  X  I N  ,  MN  »  YQ ) 

ALV{  I.  J)=YOUT(  J)/12.0 

U(  I  ,  J)=S0RT(2.  0=5'GAMA'^SIN(THETAVd)  )*ALV(  I,  JJ/RHO) 

3  CONTINUE 

Qd)  =0.5^RH0=i^U  (  I,MY)*=i^2 
C 

C      CALCULATE  DISPLACEMENT  THICKNESS 
C 

DO  31  J=1,MY 

URATIO( I  ,J)-U(  I , J)/U(  I,MY) 

F( J)=1.0-URATIO{ I, J) 
31    CONTINUE 

CALL  INTEGIF, MY, Y, RESULT] 

DELSTR{ I )=RESULT 

DELRAT(  I  )  =  DELSTRd  )  /DELSTRC  1) 
C 

C      CALCULATE  MOMENTUM  THICKNESS 
C 

DO  4  J=1,MY 

F(  J)=URATIO(  I  ,  J)*(1.0-URATI0d  tJ)  ) 

4  CONTINUE 

CALL  INTEG(F,MY,Y,RESULT) 

THETA( I )=RESULT 

THRATC I )=THETA( I)/THETA( 1) 
C 
C      CALCULATE  SHAPE  FACTOR  IN  MOMENTUM  INTEGRAL 


SHAPE( I)=DELSTR( I )/THETA( I ) 
C 
C      INTERPOLATE  P-STATIC  MANOMETER  READINGS 


DO    5    J=ltMY 

X=XOUT{ I  ,J) 

YOUT (J )=PIF2{K»XIN,MN,YPS) 

ALPd  ,  J)  =YOUT(  J)  /12.a 

P{  I, J)  =  PATM( I )-  GAMA*SIN(THETAP( I ) )*ALP{ I  ,  J) 

5  CONTINUE 

6  CONTINUE 
C 

C  CALCULATE    RESIDUAL     IN    MOMENTUM    INTEGRAL    EQUATION 


C 


DO    7    I=1,MX1 

DELX=XPOS( I+1)-XP0S (I ) 

MY=NUM( I ) 

IF(NUM( I ).GT.NUM{I+1) )       CALL     FI LL ( NUM , I t U ) 

El  =  (THETAd  +  l)-THETA(  I  )  )/DELX 

E2=THETA(I+1)^( SHAPE ( I+lJ+2.0) 

E3=(U( 1+1,MY)-U( I,MY))/(DELX*U{ I+1,MY) ) 

EPSILd+l)  =E1  +  E2*E3 

CONTINUE 

EPSILd  )=0.0 


58 


C      CALCULATE  LOCAL  TURBULENT  SHEAR  COEFFICIENT  IN  WAKE 
C 

DO  10  1=1, MXl 

DELX  =  XPOS( I  +  l) -XPOS  (I ) 

MY=NUM( I+l) 

IF(NUM( I ).GT.NUM(I+1) )   CALL  FI LL { NUM , I , ALV) 

DO  10  K=1,MY 

DO  9  J=1,K 

A( I+lt J )=2.0*JRATI0( I+1,J)-URATI0( I+1,K) 

B(I  +  1.J)=(URATI0(I+1,J  )-URATID(  I, J) )/DELX 

C1=URATI0( I+l, J) /U( I+l ,MY) 

C2=(U( I+1,MY)-U( I,MY))/DELX 

C(I+1, J)=C1^C2 

D1=SIN(THETAP(  I  +  l) ) /SI N ( THETAV ( I  +  l)  ) 

D2  =  (  ALP(  I,  J}-ALP(I  +  1,J  )  )/  (DELX=^ALV(  I  +  1,MY)  ) 

D{ I+l, J)=D1^D2 

F{  J)=A(  I+1,J)=<'(B(I  +  1,J)+C(  I  +  1,J))+D(  I  +  l,  J)/2.0 
9  CONTINUE 

CALL  INTEG(F,K,Y, RESULT) 

CTAU{ I+1,K)=2.0^RESULT 
10  CONTINUE 

WRITE(6,79) 

DO  14  1=1, MX 

WRITE(6.80)   I, MX 

WRITE(6,81)     DELSTR( I ) ,THETA{ I ) , SHAPE { I) tEPSILCI ) 

WRITE(6,82  ) 

MY=NUM( I ) 

DO    13    J=1,MY 

WRITE{6,83)    J,U(I,J  ),URATIO( I, J  )  ,P(  I,J)  ,CTAU(  I, J) 

13  CONTINUE 

14  CONTINUE 
WRITE(6,84) 
WRITE(6, 85) 
DO    15     1=1, MX 

XD  =  12.Q=^XP0S(  D/DIAMTR 

WRITE  (6,  86)  I,  XD.DELSTRd  )  ,DE  LR  AT(  I  )  ,  THET  A  ( I  )  ,  THR  AT  (  I  ) 
1,SHAPE(  I), EPSIL (  I  ) 

15  CONTINUE 
C 

C  CALCULATE    EDDY    VISCOSITY    AND    REYNOLD'S    MIXING    LENGTH 

C 

WRITE(6,87) 

DO    17    I =2, MX 

WRITE(6,88)     I 

MIN=( {NUM{ I )-l )/2)+10 

MAX=NUM(I )-10 

DO    16    J=MIN,MAX 

E1=ABS{CTAU(  I, J)*Q(  I) ) 

E2  =  A3S(U(I  , J  +  1 )-Ul I  ,J-1 )) 

E3=2.0-H 

EDDYt I, J)=( E1/RH0)/(E2/E3) 

AMIXL(  I  ,J)  =SQRT(E1/(RH0*(E2/E3)*'!^2)} 

WRITE (6,  89)   J,EDDY(  I, J) ,AMIXL(I ,J} 

16  CONTINUE 

17  CONTINUE 
C 

C      PLOT  PACKAGE  BEGINS  — MAIN  PROGRAM  ENDS. 
C 

AA  =  0. 
BB=0. 

IF     (KK.EQ.3)     AA=2.0E00 
IF    (KK.E0.2)        BB=1.0E01 
DO    100     1=1,  MX 
100       XPOSd  )  =12.-XP0S(I  ) /DI  AMTR 

CALL    DRAW (MX, XPOS, SHAP E , 1 , 5 . L AB EL3 ,T ITLE 1 , 0. , AA, 0 , 0, 0, 
13, 9, 6,0, LAST) 

CALL    DRAW(MX,XP0S,THRAT,2,4,LABEL2,TITLE2  ,0,0,0,0,0, 
10, 9, 6,0, LAST) 

CALL    DRAW( MX ,X POS , DELR AT ,3 ,3,LABEL1 ,TITLE1, 0,0, 0,0,0, 
10, 9, 6,0, LA ST) 
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WRITE(6,91)       LAST 

DO    102    I=1tMX 

MY=NUM( I ) 

IF(I.EO.l)     K=l 

IF(  I.GT.l.AND.  I  .LT..MX)     K=2 

IF(I.EO.MX)    K=3 

DO    101    J=1,MY 

UR(J)=U( I, J) 

AXIS(J)=XOUT{I ,J )/DIAMTR 

101  CONTINUE 
LABEL4=LABEL(  I  ) 

CALL     DRAW(MY,AXIS,UR,K,0,LABEL4tTITLE4tO. ,BB,0,0,OtO, 
19,15,0,LAST) 

102  CONTINUE 
WRITE{6,91)       LAST 

1000  CONTINUE 
STOP 
END 

SUBROUTINE  I NTEG { F , K,DELX, RESULT) 
IMPLICIT  REAL'^^  (A-H.O-Z) 
DIMENSION  FCl) 
;      INTEGRATION  BY  SIMPSONS  RULE 
K1  =  K'-1 
A=F(1)+F(K) 
B=0.0 

DO  1  1-2, Kl, 2 
SUM=4.Q^F( I) 
B=B+SUM 

1  CONTINUE 
C=0.0 

DO  2  1=3, K, 2 
SUM=2.0*F( I) 
C=C+SUM 

2  CONTINUE 
RESULT=DELX^(A+B+C)/3.0 
RETURN 

END 

FUNCTION  PIF2 (X,XLIST,N,FLIST) 
DIMENSION  XLISTC 1) ,FLI ST{ 1) 
BLIF(P,0,R,S,T )=( ( 0-P)*{S-T)/(R-Q)+S) 
IFIX-XLIST(N)  )  2,2tl 

1  I=N-1 
GOTO    5 

2  IF( X-XLIST( 1))     4,4,6 

4  1  =  1 

5  K  =  l 
GOTO    30 

6  K=2 

7  DO    8     1=1 ,N 

IF( X-XLIST(I)}     9,9,8 

8  CONTINUE 
I=N 

9  1=1-1 

30  BLIF1  =  BL  IF(X,XLIST{ I) ,XLIST( I  +  l) ,FLIST(  I )  ,FLIST( I  +  l) ) 

10  IF(K-l)  11,11,12 

11  PIF2=BLIF1 
RETURN 

12  IF((I+2-N))  13,13,16 

13  IF((1-1J-1)  15,14,14 

14  IF{A3S(XLIST( I-l  )-X)  -  ABS { XL  I  ST (  1  +  2 )-X ) )  16,15,15 

15  L=I+2 
GOTO  17 

16  L=I-1 

17  BLIF2  =  BLIF(X,XLIST(I)  ,XLIST{L)  ,FLIST(n,FLIST(L)J 
PIF2=BLIF(X,XLIST(I+1),XLIST{L),BLIF1,BLIF2) 

18  RETURN 
END 
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SUBROUTINE  DEC  I DE ( X ,M» N ) 
COMMON/A/  H 
DIMENSION  X(I) 
10    FORMATC///'  '  tTlO, ' NUMBER  OF  POINTS  IN  WAKE=«,I4) 
A=I.O/H 

N=IFIX( (X{M)-X(1 ))*A) 
L  =  l 

DO  5  K=ltN 
GOTO  {lt2),  L 

1  L  =  2 

5     CONTINUE 

C      N  IS  AN  ODD  NUMBER OK 

WRITE(6.10)   N 
RETURN 

2  L  =  l 

IF(K-N)  5,3,3 

3  N=N+I 

C      N  WAS  EVEN  AND  WE  FORCE  IT  TO  BE  ODD 

X(M)=X( M)+H 
C      THE  EXPERIMENTAL  DATA  SHOULD  BE  NEARLY  THE  SAME  FOR 
C      THIS  POINT  AS  FOR  THE  ORIGINAL. 

WRITE(6,10)   N 

RETURN 

END 

SUBROUTINE    F IL L ( NUM , I , U ) 

IMPLICIT    REAL*4     {A-H,0-Z) 

DIMENSION    NUM( 1) ,U(12,1} 

M=NUM( I ) 

N=NUM(I+1) 

DO     1    J=N,M 

U(  I  +  l, J )=U{  I+l ,N) 

1  CONTINUE 
WRITE(6,2)        M,N 
WRITE{6,3)        (J  (  I+l, J) , J  =  N,M} 

2  FORMATC 'J'  ,T2J , ' M='  ,13  ,1JX,' N=«  ,13//) 

3  FORMAT ( '0' , I2G10.4/) 
RETURN 

END 
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